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Abstract 

We consider a Schrodinger operator with a non-vanishing radial mag- 
netic field B = dA and Dirichlet boundary conditions on the unit disk. We 
assume growth conditions on B near the boundary which guarantee in particular 
the compactness of the resolvent of this operator. Under some assumptions on an 
additional radial potential V the operator H% — V has a discrete negative spec- 
trum and we obtain an upper bound of the number of negative eigenvalues. As a 
consequence we get an upper bound of the number of eigenvalues of H% smaller 
than any positive value A, which involves the minimum of B and the square of 
the L 2 -norm of A(r)/r, where A(r) is the specific magnetic potential defined as 
the flux of the magnetic field through the disk of radius r centered in the origin. 



1 Introduction 

Let us consider a particle in a domain Q in M 2 in the presence of a magnetic field B. 
We define the 2-dimensional magnetic Laplacian associated to this particle as follows: 
Let A be a magnetic potential associated to B ; it means that A is a smooth real one- 
form on Q C M 2 , given by A = Ylj=i a jdxj, and that the magnetic field B is the 
two-form B = dA. We have B(x) = h(x)dxi A dx 2 with b(x) = <9ia 2 (x) — <9 2 ai(x) . 
The magnetic connection V = (Vj) is the differential operator defined by 

9 

V ^ = dx-- mj - 
The 2-dimensional magnetic Schrodinger operator Ha is defined by 

3=1 



*Institut Fourier, francoise.truc@ujf-grenoble.fr Unite mixte de recherche CNRS-UJF 5582, BP 74, 
38402-Saint Martin d'Heres Cedex (France) 



1 



The magnetic Dirichlet integral Iia = {Ha-\) is given, for u G C£°(Q), by 

r 2 

h A (u) = / 2^|Vj-M| 2 |da;| . (1.1) 



3=1 



From the previous definitions and the fact that the formal adjoint of Vj is — Vj, it is 
clear that the operator Ha is symmetric on Cq°(Q). 

In [[51 we discuss the essential self-adjointness of this operator. The result in dimension 
2 is the following 

Theorem 1.1 Assume that dQ is compact and that B{x) satisfies near dVL 

b(x) > (D(x))- 2 , (1.2) 

then the Schrddinger operator Ha is essentially self-adjoint. (D(x) denotes the dis- 
tance to the boundary). This still holds true for any gauge A' such that dA' = dA = B. 

We have, using Cauchy-Schwarz inequality, 

\(b(x)u,u)\ = |([Vi, V 2 ]u,u)\ < ||Vi«|| 2 + ||V 2M || 2 u G C™(Q). 
This gives the well-known lower bound 



Wu g C °°(fi), h A (u) > 



h(x)\u\ 2 \dx\ 



(1.3) 



In this paper, we do not use the conditions (11.21) but we assume nevertheless that b(x) 
grows to infinity as x approaches the boundary. The operator H% defined by Friedrichs 
extension of the quadratic form Ha has a compact resolvent. By analogy with magnetic 
bottles on the whole space (see ID|4l[T9l), such an operator is called a magnetic bottle 
on the disk. 

We will deal with spectral estimates for the operator H%, using a perturbative 
method: introducing an additional non-negative bounded and radial potential V, we 
obtain an upper bound of the number N(A, V) of negative eigenvalues of the operator 
Hj{ — V (Theorem 12.11 ) and deduce, for any A > 0, an upper bound of the number 
N(H®, A) of eigenvalues of the operator H® smaller than A (Theorem 12.21) . Theo- 
rem [2J] can be seen as a magnetic version of the Cwikel-Lieb-Rosenblum inequality 
(see BH [TH [HI). The CLR inequality provides a bound on the number of negative 
eigenvalues of Schrodinger operators in ~§t d for d > 3 (without magnetic field) and is a 
particular case of Lieb-Thirring inequalities (see [U~5l [TTI ). 

Eigenvalue bounds were recently studied for magnetic Hamiltonians on M 2 , for 
constant magnetic fields (see IfTOl ), for Aharonov-Bohm magnetic fields (see [|3l [T4l) 
and for a large class of magnetic fields (see [[[2)). However, in [[T2l . the total magnetic 
flux 0=2^ Jr2 b(x)dx has to be finite and the dependence on the magnetic field is 
not explicit even in the radial case. In our result, the total flux is not necessarily finite 
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(see example 12.41) and the upper bound involves explicitly the square of the magnetic 
potential. 

Magnetic Lieb-Thirring inequalities were also obtained for Pauli operators (see [0 
[H), and links between magnetic and non-magnetic Lieb-Thirring inequalities were 
discussed in [0. 
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2 Main results 

We consider a smooth magnetic field B = b(x)dxi A dx 2 and a scalar potential V on 
the unit disk f2 = {x = (xi, x 2 ) G R 2 | x\ + x\ = r 2 < 1} so that 

• (Hi) K = mf xen b(x) > and b(x) -> +oo as D(x) -> 
(i.e as x approaches the boundary.) 

• (H 2 ) B is radially symmetric ( consequently we write b(r) instead of h(x)) 

• (H 3 ) V G L 1 ^), V radial and non-negative, V bounded from above . 

From assumption (Hi) and from inequality (11.31 ) we deduce that for any gauge A 
associated to B, the operator H® has a compact resolvent, and assumption (H 3 ) entails 
that the negative spectrum of H® — V is discrete, where H® — V denotes the operator 
defined by Friedrichs extension of the quadratic form Ha — V. 

Using assumption (H 2 ) we introduce polar coordinates (r, 9), (r G IR + ,# G [0,27r[) 
and consider the following magnetic potential : 

i r 

A = —a(r) sin Odxi + a(r) cos 9dx 2 , a(r) = - b(t)tdt . (2.1) 

r Jo 

We have dA = B and 

rr 

A = A(r)d9 with A(r) = ra(r) = / b(t)tdt . (2.2) 

Jo 

A(r) is the flux of the magnetic field through the disk of radius r centered in the ori- 
gin. The function a(r) = A(r)/r is well-defined (and smooth) at the origin and it is 
the amplitude of the magnetic potential A in cartesian coordinates . 
The first theorem provides an upper bound of the number N(A, V) of negative eigen- 
values of the operator H% — V where A is the magnetic potential defined by ( 12.2b . 
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From now on, A denotes this specific potential. 

Noticing that we have N(A', V) = N(A, V) for any gauge A' so that dA' = dA = B, 
we will prove the following 



Theorem 2.1 If assumptions (i7i)(f/" 2 )(if 3 ) are verified and if moreover 

3 



b(x) < M(D(x))- p , < < 



(2.3) 



for some M > 0, then 



N(A,V) < 



1 



[(- 



\-V(r)]rdr+2 / 1 + | log[r Vif ] | V(r)rdr 



for any a G]0, 1[. 

This inequality still holds when we replace in the left-hand side N(A, V) by N(A', V), 
where A' is any gauge verifying dA' = dA = B. 

The second theorem is a consequence of the first one and provides an explicit upper 
bound of the number N(H®, A) of the eigenvalues of H% smaller than any positive 
value A : 

Theorem 2.2 If assumptions (Hi) and (H 2 ) are verified and if moreover 

b(x) < M(D(x)Y 13 , < < - 

for some M > 0, then the number of eigenvalues of the operator smaller than A 
satisfies, for any a e]0, 1[, the following inequality 



N{H%,X)<c K \ 



A 



+ 



VT^ f 1 (A(r) 



a 



with 



c K 



ck 



3-1okK 



1 + log K 1 

1 

2 K 



2V1^ 

if < K < 1 
if K > I, 



rdr 



(2.4) 



Inequality A2.4\) still holds when we replace in the left-hand side N(H%, A) by N(H®, A), 
where A ; is any gauge verifying dA' = dA = B. 



Remark 2.3 The minimum of the right-hand side is obtained by taking 

a x ■ 

2 



-3/ + VP + 4/A 
A-2J 



with I :-- 



1 MM 



rdr. 



Example 2.4 Consider a magnetic field B as in the definition < \3.2\) below, and assume 
b(r) = 1 and (3 = 1. Then ck = §, the chosen gauge is A(r) = f Q b(t)tdt = 
— ln(l — r) — r and the corresponding value of I is 

I - f H1 - r)+r 'V = 2 C (3) - | . ,2.5) 
Jo r 1 

3 Proofs 



3.1 Proof of Theorem 2.1 



Let us introduce the polar coordinates x = (r,9),r G R + ,6> G [0, 2ir[. We have 
denoted by A the following vector potential : 

A = A(r)d6 with A(r) = ra(r) = / h(t)tdt . (3.1) 

Jo 

Due to assumption ( 12.31) the magnetic field we consider is of the type 

b(r) 3 
b(r) = , , a , with max6(r)<M and 5 < -. (3.2) 

v ; (1 -r)P ' [o,i[ v ; ~ P 2 v 

We first prove the following 

Lemma 3.1 If B satisfies A3.2\) . then we can find some constant C so that A writes 
A = A{r)d6 = ra(r)d8 where 

• if B < 1 max a(r) < C. 

[0,1[ 

• if = 1 a(r) = a(r) ln(l — r), with maxa(r) < C . 

Qj\T) 

• if (3 > 1 a(r) = -r— -, with max a(r) < C. 

(1 — r)P~ l [o,i[ 

7n particular j ( ) rc/r < oo . 

Proofi- 



o 



Let us explain the case > 1. The method for the case /3 = 1 is the same. 
From (13.21) we get 



< - / b(t)tdt < - / b(t)t(l - ty^dt < Ml (1 - ty^dt < 
r Jo r Jo Jo 



'1 _ r ) -|8+1 
M- — and the result follows. 

P 1 

The case f3 < 1 is straightforward. 
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□ 

We come now to the proof of Theorem l2.ll following the method of |[T3~1 . The quadratic 
form associated to H A — V can be rewritten as 



h A y(u) 



1 p2n 

Jo 



du 
dr 



V(r)\u 2 \ + r" 



[- 



iA(r)]u 



rdrdO (3.3) 



for any u G C£°([0, l[x [0, 27r[). Changing variables r = e* and denoting iu(t, 
u(e*, 9) for t g] — oo, 0[ and 9 G [0, 27r[ we transfer the form h A y(u) to 



p2n 
■oo Jo 



.dW.y ~. . . n. 



[| -*(«)]« 



(3.4) 



with 



V{t) = e*V{e% fit) = A(e*) . 



By expanding a given function iu G C^°([— oo, 0[x [0, 27r[) into a Fourier series we 
obtain that h A y(w) = ®iezh A v (we) with 



h A y{v) 



dv 
'dt' 



(£-/(t)) 2 -\/(t) 



|u 2 | tit , 



and = ILj(u>) where IL; is the projector acting as 

U e {w)(r,9) = — ^ e il ^-^w(r,9')de' . 
2tt Jo 

We write, for any a g]0, 1[ and any leZ* 



, dv l0 
: ' M ' + 



- -)/ 2 W - V(t) + (1 - a)< 2 



\v 2 \ dt . 



Let us denote by L a the operator associated via Friedrichs extension to the quadratic 
form 



Qa(V) 



dv 
'<9t ' 



(i - -)f\t) - V{t) 

a 



\v 2 \ dt . 



L a and q a depend on V and A but we skip the reference to V and A in notations for 
the sake of simplicity. Since 

h Ay > q a + (1 - a)f , 

the number N{h AV ) of negative eigenvalues of h e AV is smaller than the number of 
negative eigenvalues of L a + (1 — a)£ 2 . Hence denoting by {— /x£ } the negative eigen- 
values of L a and by ^ the set {k G N; -n% + (1 - a)£ 2 < 0} for any I G Z*, we 
get 

^(a v) < £ E 1 + N ( h Xv) ■ 
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Noticing that the sum in the right-hand side is taken over the (£, k) so that < \£\ < 
-it — we write 



N(A, V) < 



\J\ — a 



(3.5) 



fcSN 



Let us extend the functions / and V to R by zero and denote respectively by fi and V\ 
these extensions. 

Since C£° ([— oo,0[) C Cq°(R), the negative eigenvalues {—v%} of the operator L° 
associated via Friedrichs extension to the quadratic form 



dv 

'at' 



(i - -)A 2 (0 - 

a 



\v 2 \ dt 



verify 



Ev^<E^ (3.6) 

fceN fceN 

Applying the sharp inequality of Hundertmarkt-Lieb-Thomas ifTTI (see Appendix ) to 
the operator we get 



E 

feeN 



< 



;--i)^( f ) + K(t) 



dt 



o r 



(--1)/ 2 W + ^) dt 

a 



< 



a 



r<ir . 



To conclude we need the following 

Lemma 3.2 Assume that K = inf^n h(x) > O.Thenfor any e e]0, 1[ 



N(h% v ) = N(h% -V)< 



K l-e)K .. 
1 + | lo S (\/ —r)\ 



In particular 



i r 



N(h° AV )<2 l + \\og(VKr)\ V(r)rdr . 



(3.7) 



V(r)rdr . (3.8) 



(3.9) 



Proof.- 



1 



Step 1 : From (OT) we get that h A (u) > K J n \u\ 2 \dx\ Mu e 
Co°(f2), which implies for h A (returning to the variable r and con- 
sidering V = 0), 



<9u> 
<9r 



rdr 



> A" f \w\ 2 rdr Vu> G C °°([0, 1]) . 
Jo 

We write for any e e]0, 1[ 

JV(/i^ -V) < JVft+Mif-V) < iV (h% + - ^ 

(3.10) 

where we have used the fact that multiplying an operator by a posi- 
tive constant does not change the number of its negative eigenvalues. 

Step 2 : We establish the following upper bound : 

N(h° A0 + l-V) = N(h% v + l)< [ [1 + | log r\]V(r)rdr. 

Jo 

(3.11) 



We have 



h%y( W ) 



I ^ 

dr 



rdr 



> 



o 



— | -V{r)\w 



+ [r- 2 A 2 (r) - V(r)] \w 

rdr G Cg°([0, 1]) . 



By the variational principle, 

N(h% v + 1)<N(P + 1-V), (3.12) 

where P is the operator generated by the closure, in L 2 ([0, 1], rdr) 
of the quadratic form 

, dw , 



o dr 



'rdr, w e C °°([0, 1]) 



Considering the mapping U : L 2 ([0, 1], rdr) — > L 2 ([0, 1], dr) de- 
fined by (Uf)(r) = r 1 / 2 /(r) we get that 



iV(P + 1 - V) < N{T + 1-V) 



(3.13) 



where the operator T = UP U 1 is the Sturm-Liouville operator on 
L 2 ([0, 1], dr) acting on its domain by 



(T oM )(r) = - M "(r) 



it(r) 
4r 2 
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, n(0) = n(l) = 0. 



(3.14) 



The upper bound (13.1 II ) will follow from the properties of G(r, r, 1), 
the diagonal element of the integral kernel of (T + Precisely 
we have 

G(r,r, 1) < r(l + | logr|), r 6 [0, 1[ . (3.15) 

The proof of (13.151 ) is given in Appendix . The Birman-Schwinger 
principle then yields 

N(T + 1-V)< [ G(r,r,l)V(r)dr < [ [1 + | log r|] V(r)rdr . 
Jo Jo 

(3.16) 

This ends the proof of (13.111) . together with the inequalities (13.121) 
and (I3T3T) . 

Step 3 : We mimick the previous method to get, for any strictly pos- 
itive number k 

N(h° A0 + k 2 - V) < [ [l + \\og(kr)\]V(r)rdr . (3.17) 
Jo 

Due to the Birman-Schwinger principle it suffices to prove that, for 
any strictly positive number k 



G(r,r,k 2 ) < r(l + |log(»|), r e [0,1[. 
This is done in Appendix . 



(3.18) 



Step 4 : Returning to (13.101) and applying (13.171) with k 
and — instead of V we get, for any e e]0, 1[ 



_ (l-e)K 

e 



N(h% - V) < N (h° A 



1 - e)K _ V 

e e 



< 



1 + I 10E 



(l-e)K 



r)\ 



and taking e = \ we obtain Lemma [3T2l 



(3.19) 

V(r)rdr , 
(3.20) 



□ 

Theorem 12. 1 1 follows from Lemma \32\ together with inequalities (13.51) . (13.61) . and 
(EZ]). 



3.2 Proof of Theorem 2.2 



Noticing that for any A > the constant potential V(x) = A is in L 1 (f2), and that 
N(A, A) denotes the number of eigenvalues of the operator less than A, we apply 
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Theorem l2.1l to V(x) = A. To get the result it suffices to compute J Q [1 + | log(fcr) |] rdr 
We get after computation that 



i 

[1 + | log(Ax) |] rdr = 7fc , (3.21) 



with 



3-21ogfc 
Ik = t if k < 1 

l + 21ogfc 1 

7fc = 2 h — T77 if fc > 1 

' 4 2fc 2 



3.3 Proof of Remark 2.3 



To get the minimum over the values of a we study the sign of the expression, for any 
a e]0, 1[, of 

A VI - " r 

^a(«) := ~ 7= + / • 

2V1 — ol a 

A direct computation shows that the value a\ which realizes the minimum of g\(a) is 
the positive solution of 

a 2 (A-2/) + 6a/-4/ = 0. (3.22) 



4 An asymptotic eigenvalue upper bound 

From Theorem I2.2l we get easily an asymptotic estimate for the right-hand side of (12.41 ) 
when A tends to oo : 

Corollary 4.1 If assumptions (Hi) and (H 2 ) are satisfied and if moreover 

h(x) < M(D{x)Y P , < < | 

for some M > 0, then the number of eigenvalues of the operator smaller than A 
satisfies, as A — > oo 

N(H%, A) < + c K )X + + 0(1) , (4.1) 



where 

I= f 1 f A (ry 

Jo 

and 



rdr 



3 -log ^ ., _ ^ u . . 

c K = if < A < 1 
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c K 



1 + log K 1 
2 + K 



if K>1 



Inequality < \4.1\) still holds when we replace in the left-hand side N(H^, A) by 
N(H®, A), where A' is any gauge verifying dA' = dA = B. 

Proof. - 

We define as previously, for any a e]0, 1[, 



A y/l-a 
2y/l - a a 

and we want to determine the asymptotic behavior as A tends to oo of 

g\(a\), where a x is the minimum of g\{a). 

From (13.221) we compute the following asymptotic s 

v / r^:=i-^ + o(-), 



and this gives the result. 



□ 



Remark 4.2 The leading term in the estimate d4.il) is of the same order than the lead- 
ing term in the Weyl formula for the Dirichlet Laplacian (corresponding to the case 
A = ) in the unit disk. 

5 Appendix 

5.1 The inequality of Hundertmarkt-Lieb-Thomas 

We recall the sharp inequality of Hundertmarkt-Lieb-Thomas [fTTTl 
Theorem 5.1 Let 

Lv(t) = -v" (t) - W(t)v(t), W>0 W e L X (R) 

be defined in the sense of quadratic forms on R, and assume that the negative spectrum 
of L is discrete. Denote by {—Uk,k G N} the negative eigenvalues of L. Then 



1 r+°° 

J2v^<- w(t)dt. 
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5.2 The Green function G(r, r', 1) of the operator T . 

Let us compute the diagonal element for the Green function G(r, r', 1) of the operator 
T defined by (13.141) . G(r, r', 1) is the solution of 

((T + l)u)(r) = MO, u(0)=m(1) = 0. (5.1) 

We have 

G(r, r', 1) = A\U\{r) + A 2 u 2 (r) r < r' 

G(r, r', 1) = B\U\{r) + B 2 u 2 {r) r > r' , 
where u±(r) = y/rl (r) andu 2 (r) = y/rKo(r) are independent solutions of the related 
homogeneous equation, (7 and K are the modified Bessel functions). 
The coefficients depend of r' but we omit the indices for the sake of clarity. Due to the 
boundary conditions and to the fact that the derivative (with respect to r) of G(r, r', 1) 
has the discontinuity in r' of a Heaviside function, they satisfy : 

Aiui(0) + A 2 u 2 (0) = Siui(l) + B 2 u 2 (l) = 

Ft A 'Mr') R A _ Mr') 

A l — TTT, A tf 2 — A 2 — 

W(r') W(r') 
where W(r') is the value of the Wronskian of U\ and u 2 taken at the point r'. 
The first equation is always satisfied since wi(0) = u 2 (0) = 0. Let us set A 2 = 0. We 
have W(r') = M , 1 (r')M 2 (r / ) — ui{r')u' 2 (r') = r'W(r') where W(r') is the Wronskian 
of the modified Bessel functions J and K . As r'W(r') = 1 (see [2]), we get after 
solving the above system, and doing r = r' : 



G(r, r, 1) = U\{r) 



u 2 (l) 

~ u i{ r ) 7TT + u 2 [r) 
uAl) 



= rl {r) 

Using again the properties of the modified Bessel functions (see (El) we can write 

G(r,r,l)<rI (r)K (r). 

The function 

I (r)K (r) 

g{r) 



1 + I logr| 

has a limit at r = equal to 1 (see 0), so 

h(r)K (r) 

c = max — — f (5.2) 

[o,i[ 1 + | log r | 

exists and 

G(r,r, 1) < cqt{1 + | log r|), r e [0, 1[. 

Numerics suggest that g is decreasing on [0,1], so that one should have c = 1. In next 
subsection, we give the proof of this result, which can not be found to our knowledge 
in the literature, and has been communicated to the author by J. P. True [|20l : 
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Proposition 5.2 Vr e]0, 1] : Io ^ K °^ < 1. 

1 — log r 



5.3 Proof of Proposition 1572 



The modified Bessel function I can be written as 

W = E^ = 1 + I + - (5-3) 

Therefore we have 



+ °° (lf.) fc 2 

1 < h{r) <^H^ = e~ 



k=0 

and 

Vr e]0, 1] : 1 < J„(r) < . (5.4) 
According to the expression of the modified Bessel function K 

Mr) = -(log(r/2) + 7 )/o(r) + E (E < 5 - 5 ) 

fe=i i=i 3 

where 7 denotes the Euler constant, we compute that 

K (r)I (r) - (1 - logr) = S(r) - 1 , (5.6) 
where S(r) denotes the following function : 

+00 k K (r 2 \k 

^ ( ^ iXj 

\^ j) k\ 2 

k=i 3=1 J 

Proposition l5.2l is then a straightforward consequence of the following Lemma 
Lemma 5.3 



5(r) = (1 - J (r) 2 ) logr - ( - log 2 + 7 )/ (r) 2 + J (r) E ( E ~) " < 5 - 7 > 



Vr e]0,l] : 5(r) < 1. 

Proof. - 

The function 5(r) splits into 3 positive parts, which we study separately . 
• An upper bound for (1 — Io(r) 2 ) logr. 

2 

From (15.41 ) we deduce 1 — h(r) 2 > 1 — e^, and : 
Vre]0,l] : < (1 -I (r) 2 ) logr < (e^ - logr) < 0,11. 
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An upper bound for ( — log 2 + 7 j Io(r) 2 . 

A straightforward computation gives — 7 + log 2 < 0.12 so using that 

-fo( r ) < eI we g et 

( -log 2 + 7) J (r) 2 < 0.16. 

+00 ^ k ^ j-rf^fc 



• An upper bound for J (r) ^ ^ ^ . 



fc=i j 

k 



For fc G N*, we set s fc = -. We have si = 1 . For k > 2, 



according to the inequality 



1 f k dt 
< / — = log k — log (A; — 1). 



we get that: 

k 



^2 - < log A; 



and for any integer fc , s fe < 1 + log k. Thus 

v r f ixi! < v r 1+logfc > i ( ^ )fc 

^ 7/ A;! 2 - ^ V fc! / A;! 

k=l j=l J k=l 

Noticing that, for any integer k > 1 

1 + log k 1 + log fc 
" A;! " fc 



we can write, Vr e]0, 1] : 



+°° , k ^ (rt\k +°° (rt\k 2 



fe=i 3=1 J k=i 
Finally we have, for any r g]0, 1] 

+00 k 1 / t 2 \k 

W E (E 7) %fSei( e J- 1)^ 0.364. 

fc=l J=l J 

Summing the 3 previous estimates one gets : Vr e]0, 1] : 8{r) < 0.11 
0.16 + 0.37 < 1 . 



The optimality of the value c = 1 is due to the fact that 

lim g.(r)J.(r) = L 
r->o+ 1 — In r 



□ 
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5.4 The Green function G(r, r', k 2 ) of the operator T 

We now compute the diagonal element for the Green function G(r, r',k 2 ) of the oper- 
ator T defined by (13.141) . G(r, r', k 2 ) is the solution of 

((T + k 2 )u) (r) = 8 r ,(r), u(0) = u(l) = . (5.8) 

We have, as previously 

G(r, r, k 2 ) = «i(r) 

where Ui(r) = y/rIo(kr) and U2(r) = y/rK (kr) are independent solutions of the 
related homogeneous equation. This leads to 

G(r, r, k 2 ) = rl (kr) 
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